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Abstract 

The B c — > B q P, B q V decays are studied with the QCD factorization approach (where P and V 
denote pseudoscalar and vector mesons, respectively; q = u, d and s). Considering the contributions 
of both current-current and penguin operators, the amplitudes of branching ratios are estimated 
at the leading approximation. We find that the contributions of the penguin operators are very 
small due to the serious suppression by the CKM elements. The most promising decay modes are 
B c — ► Bs vr, B s p, which might be easily detected at hadron colliders. 
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I. INTRODUCTION 



The B c meson is one of the unique "double heavy-flavored" binding system in the standard 
model (SM). The study of the B c meson has received a great interest, due to its special 
properties: (1) The B c meson carries open flavors. We can study the two heavy flavors of 
both b and c quarks simultaneously with the B c meson. (2) The B c meson can serve as a great 
laboratory for potential models, QCD sum rules, Heavy Quark Effective Theory (HQET), 
lattice QCD, etc. (3) The B c meson has rich decay channels, because of its sufficiently large 
mass and that the b and c quarks can decay individually. The B c meson decays may provide 
windows for testing the predictions of the SM and can shed light on new physics beyond 
SM. 

The B c mesons are too massive to access at the I?-factories near T(4S). They can be pro- 
duced in significant numbers at hadron colliders. The B c meson has been firstly discovered 
by the CDF Collaboration ll| . Recently the CDF and DO Collaborations announced some 
accurate measurements 0, [J with part of their available data. Much more B c mesons and 
detailed information about their decay properties are expected at the Large Hadron Collider 
(LHC) which is scheduled to run in this year. It is estimated that one could expect around 5 
x 10 10 B c events per year at LHC % [5| due to the relatively large production cross section 
6] plus the huge luminosity C = 10 34 cm~ 2 s _1 and high center-of-mass energy ^/s = 14 TeV 
7|. There seems to exist a real possibility to study not only some B c rare decays, but also 
CP violation and polarization asymmetries. The study of the B c meson will highlight the 
advantages of B physics at hadron colliders. 

The B c meson is stable for strong interaction because it lies below the threshold of the 
B-D mesons. The electromagnetic interaction cannot transform the B c meson into other 
hadrons containing both b and c heavy quarks, because the B c meson itself is the ground 
state. The B c meson decays via weak interaction only, which can be divided into three 
classes: (1) the b quark decay (b — > c, u) with c quark as a spectator, (2) the c quark decay 
(c — > s, d) with b quark as a spectator, and (3) the weak annihilation channels. In the 
B c meson, both heavy quark can decay weakly, resulting in its much shorter lifetime than 
other 6-flavored mesons, i.e. tb c < \ T B q (where q = u, d, and s) jsj]. Rates of the Class (1) 
and (2) are competitive in magnitude. The Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements \V c b\ <C \V CS \, that is in favor of the c-quark decay greatly, whereas the phase space 
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factor m\ 3> compensates the CKM matrix elements a lot for the two flavors [loj. In 
fact, the dominant contributions to the B c lifetime comes from the c-quark decays [Class (2)] 
(~ 70%), while the 6-quark decay [Class (1)] and weak annihilation [Class (3)] are expected 
to add about 20% and 10%, respectively |4[. 

The B c meson decays have been widely studied in the literature due to some of its 
outstanding features. (1) The pure leptonic B c decays belong to the Class (3), which are 
free from strong interaction in final states and can be used to measure the decay constant 
/b c and the CKM elements \V c b\, but they are not fully reconstructed due to the missing 
neutrino. (2) The semileptonic B c decays provide an excellent laboratory to measure the 
CKM elements \Vcb\, \V u b\, \V CS \, \Vd\ and form factors for transitions of B c — > b- and c- 
flavored mesons. The first signal of B c is observed via this mode pj. The most difficult 
theoretical work at present is how to evaluate the hadronic matrix elements properly and 
accurately. (3) The nonleptonic B c decays are the most complicated due to the participation 
of the strong interaction, which complicate the extraction of parameters in SM, but they 
also provide great opportunities to study perturbative and non-perturbative QCD, final state 
interactions, etc. 

he earlier nonleptonic decays of B r meson has been studied in 
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While c-quark decays take the lion's share of the B c lifetime, the study on the Class (2) 
has not received enough attention. This can be explained by the fact that on the one hand 
the available data on the B c meson is very few, on the other hand it is assumed that the 
long distance effects and final state interferences might be quite huge and that the Class (2) 
decays were hard to detect experimentally. Accompanied by the LHC being about to run, 
the future copious data require more accurate theoretical predictions from now on. In this 
paper, we shall concentrate on the B c — > B^P, B q V (here P and V denote pseudoscalar 
and vector mesons, respectively; q = u, d and s) decays in Class (2) with QCD factorization 
approach. Now let us outline a few reasons and arguments below. 

1. From the experimental view 

• The initial and final 6-flavored mesons, i.e. the B c and B)*\ all have a long 
lifetime due to their decays via the weak interaction. Considered the relativistic 
boost kinematically due to their large momentum obtained from huge center- 
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of-mass energy, their information would be easily recorded by the multipurpose 
detectors sitting at the hadron colliders interaction regions (see pj for details). 

• Although it is perceived that the hadron collider environment is "messy" with 
high backgrounds, the B c — > B q *^P, B q V decays are measurable due to the "clean" 
final states. Since the B c meson carries charge, the final meson is tagged 
explicitly by the initial B c meson. The other light meson in the final state could 
also be identified effectively by the conservation low of both momentum and 
energy, because the dedicated detectors at LHC has excellent performance on 
trigger, time resolution, particle identification and so on (see j3] for details). 

2. From the phenomeno logical view 

• With very high statistics, we can carefully test the various theoretical models, 
precisely determine the CKM elements, and meticulously search for the signals of 
new physics. This requires more accurate theoretical predictions. In this paper, 
we shall study the B c — > Bj*^P, B q V decays with QCD factorization approach, 
including the contributions of both current-current and penguin operators. 

• In the rest frame of the B c meson, the velocity P B <*) of the B^*' meson is very 
small due to its large mass, not exceeding 0.18. The ratio of velocity /?py //? (*) 
> 5.5, which is very different from that in the two-body D meson decays where 
the ratio of velocities of final states is close to one. This may indicate that the 
final state interferences for B c — > B^*'P, B q V decays might not be so strong as 
that in D mesons. If it holds true, it will benefit us in determining the CKM 
elements V cs and V c d, the B c — > 5^ transition form factors, etc. In this paper, 
we shall neglect the effects of final state interferences for the moment. 

This paper is organized as follows: In Sec. [Til the theoretical framework is discussed. 
To estimate the amplitude of the branching ratios, the master QCD factorization (QCDF) 
formula are applied to the B c — > B q *^P, B q V decays at the leading approximation. Section 
ITTT1 is devoted to the numerical results. Finally, we summarize in Sec. [TV] 
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II. THEORETICAL FRAMEWORK 



A. The effective Hamiltonian 

Using the operator product expansion and renormalization group (RG) equation, the low 
energy effective Hamiltonian relevant to the B c — > B^P, B q V decays can be written as 

n ^ = £K 9l ^4 [CMQi + cmq 2 ] + 1 £v vq vzc i {riQ i } + H.c, (1) 

V q-^=d,s q—d,s 

q2—d,s i— 3,. ..,10 

where V uq V* q ^ is the CKM factor. The cases q = d and q = s can be treated separately and 
have the same Wilson coefficients Ci(fi). The expressions of the local operators are 



Qi = 


{u a Qia)v-A{q2f3Cp)v-A-, 


Q2 = 


(u a qip)v-A(Q2/3C a )v-A 


(2) 


Qi = 


{u a c a ) v _ A Y,%q'p)v-A, 


Q* = 


(u a Cp) V „ A J2(q'p<la)v-A, 

q' 


(3) 


Q 5 = 


q' 


Qe = 


q' 


(4) 


Qr = 


3 

-(u a c a ) V - A J2 e i'(yMv+ A i 
q' 


Qs = 


3 

~(u a Cp) V - A J2 Zq'{m' a )v+A, 

q' 


(5) 


Q9 = 


3 

-(u a c a ) v ^. A J2 e q'%Qp)v-A, 
q' 


Qw - 


3 

= ^aC^V-A^eq'iqWah-A, 

q' 


(6) 



where the summation over the repeated color indices (a and (3) is understood. The Dirac 
current (qiq2)v±A = <?i7(l±75)<?2- <?' denotes all the active quarks at scale fi = 0(m c ), i.e. 
q' = u, d, s, c. e q i denotes the electric change of the corresponding quark q' in the unit 
of |e|, which reflects the electroweak origin of Qj, ■ ■ •, Qiq. The current-current operators 
(Qi, Q2), QCD penguin operators (Q3, • • •, Qq), and electroweak penguin operators (Q7, 



• • •, Qio) form a complete basis set under QCD and QED renormalization [c 

The effective coupling constants — Wilson coefficients Ci(fi) — are calculated in pertur- 
bative theory at a high scale fi ~ mw and evolved down to a characteristic scale /i ~ m c 



using the RG equations. The Wilson coefficient functions are given by 



38 1 



C(jl) = t/ 4 (/x, ii b )M(ii b )Us{nb, Hw)C(ii w ) (7) 

Here Uf([if,fii) is the RG evolution matrix for / active flavors, which includes the RG- 
improved perturbative contribution from the initial scale fii down to the final scale /z/. The 
M(/i) is the 10x10 quark-threshold matching matrix. The corresponding formula and ex- 
pressions can be found in Ref. 



381 ] . The Wilson coefficients Ci(fi) have been evaluated to the 
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next-to-leading order (NLO). Their numerical values in the naive dimensional regularization 
(NDR) scheme are listed in Table [B 



B. Hadronic matrix elements within the QCDF framework 

For the weak decays of hadrons, the short-distance effects are well-known and can be 
calculated in perturbation theory. However, the nonperturbative long-distance effects re- 
sponsible for the hadronization from quarks to hadrons still remain obscure in several as- 
pects. But to calculate the exclusive weak decays of the B c meson, one needs to evaluate 
the hadronic matrix elements, i.e., the weak current operator sandwiched between the initial 
state of the B c meson and the concerned final states, which is the most difficult theoretical 
work at present. Phenomenologically, these hadronic matrix elements are usually parame- 
terized into the product of the decay constants and the transition form factors based on the 



argument of color transparency and the naive factorization scheme (NF) 39(. A few years 



ago, Beneke, Buchalla, Neubert, and Sachrajda suggested a QCDF formula to compute the 
hadronic matrix elements in the heavy quark limit, combining the hard scattering approach 
with power counting in I/tuq 40j (here rriQ is the mass of heavy quark). At leading order in 
the power series of heavy quark mass expansion, the hadronic matrix elements can be fac- 
torized into "non-factorizable" corrections dominated by hard gluon exchange and universal 
non-perturbative part parameterized by the form factors and meson's light cone distribution 
amplitudes. This promisin g ap proach has been applied to exclusive two-body nonleptonic 



B u , B d , B s decays 



41 



42 



43| . It is found that with appropriate parameters, most of the 
QCDF's predictions are in agreement with the present experimental data. In this paper, we 
would like to apply the QCDF approach to the B c — > B^P, B q V decays. 

In the heavy quark limit m c ^> Aq C d, up to power corrections of order of the AQ CD /m c , 
using the master QCDF formula, the hadronic matrix elements for the B c — > B^M decays 
(M = P or V) can be written as 4(]] 

(B^Mp^) = F^*' JdzH(z)$ M (z) (8) 

where F Bc ~* Bq ' is the transition form factor and &m( z ) is the distribution amplitudes for 
the meson of M, which are assumed to be nonperturbative and dominated by the soft 
contributions. The hard-scattering kernels H(z) can be calculated in the perturbative theory. 
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For details about the QCDF formula Eq.(jHJ), please refer to Ref. |40|. 

To estimate the branching ratios approximately and to have a sense of the order of 
amplitudes, we shall adopt a rough approximation, i.e. at the leading order of a s . Within 
this approximation, the hard-scattering kernel functions become very simple, H(z) = 1. That 
is to say the long-distance interactions between M and B c -B^> system could be neglected. So 
the integral of $m (2) reduces to the normalization condition for the distribution amplitudes. 
Furthermore, according to the arguments of QCDF [4^], the hard interactions with the 
spectator are power suppressed in the heavy quark limit. Therefore it is not surprisingly to 
reproduce the result of NF. 

In our paper, the annihilation amplitudes are ne glec ted due to some reasons. (1) Ac- 
cording to the power counting arguments of QCDF [40], compared with the leading order 
contributions to the hard scattering kernel, the contributions from annihilation topologies 
are power suppressed. (2) The annihilation amplitudes are suppressed by the CKM elements. 
For the decay modes concerned, the CKM factors in the non-annihilation amplitudes are 
VudV* s ~ 1, V US V* S ~ A, V u dV* d ~ A and V us V* d ~ A 2 , while the annihilation amplitudes are 
proportional to the CKM factors of V c bV* b ~ A 5 . 

The explicit expressions of decay amplitudes for B c — ► Bff'P, B q V decays are collected 
in appendix |A] In our paper, we define 

<H = Q + jj-Ci+x (i = odd) (9) 
Oj = Ci + j-j-Ci-i {i = even ) ( 10 ) 
where % runs from 1 to 10, C,- are the Wilson coefficients. N r = 3 is the color number. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

Within the QCDF approach, the decay amplitudes depend on many input parameters 
including the CKM matrix elements, decay constants, form factors, etc. These parameters 
are discussed and specified below. 
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A. The CKM matrix elements 



We will use the Wolfenstein parameterization. Phenomenologically, it is a popular ap- 
proximation of the CKM matrix in which each elements is expanded as a power series in the 
small parameter A. Up to C(A 6 ), the CKM elements can be written as [38] 



V ud 

V us 

v ub 
v cd 

v cs 
v cb 

V td 
V u 

v tb 



i_I A 2 -Ia 4 + o(a 6 

2 8 
A + C(A 6 ) 

AX 3 (p - irj) 

-A + A 2 X 5 
1 



{p + irf) +0{X 



z X 2 - -X 4 - -A 2 X 4 + 0(X 6 
2 8 2 v 



AX 2 + 0{X 



6> 



AX 



l-(p + zr / )(l--A 2 )]+0(A 



-AX 2 + AX A [--(p + ir l )\ +0{X 

1 - -A 2 X 4 + 0(A 6 ) 
2 



The global fit for the four independent Wolfenstein parameters gives 



a = 0.81818:8??, 



A = 0.2272±0.0010, 



P 



n 991 +0.064 
U-^ 1 -0.028' 



fi - n 04Q+0.017 

7 — u -°^ u -0.045 



where the relationship between (p, rj) and (p, rj) is js 

p + irj 



VI - A 2 A 4 (p + 



v/T^A^l - A 2 A 4 (p + ^)] 
If not stated otherwise, we shall use their central values for illustration. 



(11) 
(12) 

(13) 

(14) 

(15) 
(16) 
(17) 

(18) 

(19) 



(20) 



(21) 



B. Decay constants and form factors 

In principle, information about the decay constants and transition form factors of mesons 
can be obtained from experiments and/or theoretical estimations. Now we specify these 
parameters. The decay constants fp and fy corresponding to the pseudoscalar and vector 
mesons respectively, are defined by 



(^(9)1(5192)^10) = -ifpq», (V(q,e)\(q iq2 ) v _ A \0) = f v m v e*», (22) 
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where e* is the polarization vector of the vector meson V. In this paper, we assume ideal 
mixing between u and <\> mesons, i.e. cu = (uu + dd) /y/2 and (p = ss. In fact, the B c — ► 
B^'cj) decays are not possible, because the B c meson lies below the threshold of the B^'<f> 

|, adopting the 



system. As to the 77 and 7/ mesons, we take the convention in Ref. 
Feldmann-Kroll-Stech mixing scheme. Neglecting the possible compositions of both rj c = cc 
and glueball gg, the 77 and 77' are expressed as linear combinations of orthogonal states t] q 
and 77 s with the flavor structure qq = (uu + dd)/\/2 and ss, respectively, i.e. 

/ 

(23) 





where = (39.3±1.0)° [44j is the 77-77' mixing angle. So the decay constants related to the 
77 and 77' mesons can be defined by 

' fv fv \( cos< ^ - si < ^ ( 
\ fv' ft ) V sin< ^ cos< ^ 
(0|ffVM|»7 (/) (p)) = itfioP^ (0|s7 M 7 5 s|77 (/) (p)) = if° 0)Pll . 
The matrix elements of the pseudoscalar densities are defined by 




45] 



m 



2m , 



(24) 
(25) 

(26) 



The numerical values of the decay constants are collected in Table HT1 If not stated otherwise, 
we shall take their central values for illustration. 
The transition form factors are defined as 39] 



(P(k)\(q m ) V - A \B(p)) = ( p + kYF^ p {q 2 ) + 



m 



B 



iB->Pf„2 



(q 2 )'Fr P (q 2 ) (27) 



(V(k,e)\(q m ) v . A \B(p)) = t ^ qfl 2m v A^ v (q 2 )+te*(m B + m v )A^ v (q 2 ' 



C-p 



e -p 



(p + k)pAf V (q 2 ) - i^q,2m v A^ v (q 2 ^ 



m B + my 



2V 



B^V I Jl 



(q 2 



(28) 



m B + mv 

where Fq A , V and ^0,1,2,3 are the transition form factors, q = p — k. In order to cancel the 
poles at q 2 = 0, we must impose the condition 



(0) = Jf-^(O) 



aB^V 
^0 



(0) 



A B r V 



(0), 



2m v A^ v (0) = (m B + m v )Af ^ v (0) - (m B - m v )A^ v (0). 



(29) 
(30) 



In our paper, only the B c — ► B^> transition form factors appear in the amplitudes within 
the "spectator" model where the spectator is the 6-quark for the concerned processes. Their 
numerical values are collected in Table IIHI From the numbers in Table 1111} we can see clearly 
that, due to the properties of nonperturbative QCD, there are large uncertainties about the 
form factors with different theoretical treatments. Here, we notice the fact that the velocity 
of the final state B^ meson is very small in the rest frame of the initial B c meson, as that 
mentioned in Sec. [H It is commonly assumed that the velocities of the 6-quark in the rest 
frame of the 6-flavored mesons should be close to zero. The meson is neither fast nor 
small. By intuition, the overlap between the initial and final states should be huge, close to 



unity, as that argued in 37|. So for illustration and simplification, we will take the same 



value for the transition form factors, i.e. F 10 (0) = A (0) = 0.8. 
C. Quark masses 

In the decay amplitudes, there exist the "chirally enhanced" factors which are associated 
with the hadronic matrix elements of the scalar and pseudoscalar densities, for example, R c ± 
in Eq. llAip . These factors are formally of order the Aqcd/^c, power suppressed in the heavy 
quark limit, but numerically close to unity because the mass of the c quark is not infinity 
in practice. The current quark masses in the denominator appear through the equations of 
motions and are renormalization scale dependent. Their values are 



m u {2 GeV) = 3±1 MeV, m d (2 GeV) = 6.0±1.5 MeV, 

m s {2 GeV) = 103±20 MeV, m c (m c ) = 1.24±0.09 GeV. 



(31) 



Using the renormalization group equation of the running quark mass 38J , 



their corresponding values at a characteristic scale /i ~ m c can be obtained. 



D. Numerical results and discussions 

The numerical results are listed in Table HVl where Br T corresponds to the contributions of 
the current-current operators only, Br T+Ps corresponds to the contributions of both current- 
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current and QCD penguin operators, Br T+Pa+Pc corresponds to the contributions of both 
current-current and penguin operators, i.e. Qi, ■ ■ •, Qio- 

Here, we would like to point out that these numbers are just the qualitative estimations 
on the order of amplitudes, because many of the subtleties and details, such as final state 
interactions, the renormalization scale dependence, the transition form factors, the strong 
phases, and so on, all deserve the dedicated researches but are not considered here. 

From the numbers in Table IIVI we can see 

• The contributions of both QCD and electroweak penguin operators are very small 
for B c — > B(*'P, BV decays, compared with those of the current-current operators. 
This is very different from that of the B u ^ s meson decays. The reason is that the 
contributions of penguin operators are seriously suppressed by the CKM elements. 
The CKM elements corresponding to different topologies for c-quark decay in the B c 
meson are listed below. 



tree topologies 


penguin topologies 


annihilation topologies 


v ud v* 
v us v; d 


~ 1, V US V* ~ +A 
~ A 2 , V ud V* d ~ -A 


v ud v: d + V US V* ~ A 5 


VcbV* b ~ A 5 



So, for the B c — ► B^P, BV decays, the effects of new physics contributed via the 
penguin topologies might be tiny and not detectable even with large statistics, due to 
the serious suppression by the CKM elements. 

There are clear hierarchy of amplitudes of the branching ratios. According the CKM 



elements and the coefficients of a 12 


, these decay modes 


are divided into different 


listed below. 










cases 


processes 


coefficients 


the CKM elements 


order of branching ratios 


case la 


c — ■> s 


Oi 


v ud v; s ~ i 


~ io- 2 


case lb 


c —>• s 


a% 


V US V* ~ A 


~ 10^ 3 




c — > d 


ai 


v ud v: d ~ A 


~ 10" 3 


case lc 


c — > d 


a% 


v us v: d ~ a 2 


~ 10" 4 


case 2 a 


c — > u 


a 2 


v ud v; s - 1 


~ 10^ 5 


case 2b 


c — > u 


a 2 


v us v c * s , v ud v; d ~ a 


~ IO" 6 — 1(T 7 


case 2 c 


c — > u 


a 2 


v us v; d ~ a 2 


~ io- 8 
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The decay modes determined by a\ have comparatively large branching ratios, which 
should be detectable experimentally, especially the CKM favored decay modes, such 
as B c — > B^*'tt, B s p, might be the promising decay modes to be measured in hadron 
colliders. Due to the great branching ratios of the decay modes determined by ai, the 
B c mesons can be used as a source of the B s mesons if the B c is produced copiously, 



as that stated in Ref.|lOj. The decay modes determined by ai have comparatively 
small branching ratios, which are hard to detect experimentally, especially the CKM 
suppressed decay modes, such as B c — > B^K°, B U K*°, their branching ratios are too 
tiny to be measured. 

Although the B c — > B^/rf decays belong to the case 2b modes, their branching ratios 
are abnormally small, order of 10 -8 . This can be explained by the fact that on one 
hand the physical space phase available is too small, on the other hand there are 
large destructive interactions between f5a,2 and f^a-i due to the serious cancellation 
between the CKM elements V ud V* d and V US V* S . 

The relations among the B c — > B^P, B q V decay mode become very simple since the 
effects of penguin topologies is too tiny to be considered. We can use these relations 
to determine and overconstrain some parameters, such as the CKM elements, the 
form factors, etc. In addition, in estimating and measuring these parameters, the 
ratios of the branching ratios can be used to cancel and/or reduce largely theoretical 
uncertainties and experimental errors. For example 

Br(fl+-fl27r+) „ V*V« U ^ Vu d V* U _ Br(B+^B° s 7r+) 

(34) 
(35) 
(36) 
(37) 
(38) 



Br(B+- 




v us v; d f K 


V us V*f K Br{B+^ 






Br{B+- 




V*d Fo c ^ B 


1 (m 2 Br -m 2 Bd )_Br(B+ 


^B° S K + 


) 


Br{B+- 




V* F B ^ B 


'{m 2 Bc -m 2 Bs )'~Br(B+ 




) 


Br{B+- 






Br(B+^B° d K+) ^ 


f K F^- 


■B d 


Br{B+- 




fpF B^B d 


Br(B+^B° d K*+)~ 


fK*F Bc ~ 


+B d 


Br(B+- 




TpB c ^B d 


Br(B+^B° d K+) 






Br(B+- 


^B?n+)' 


.B C ^B* 


Br(B+^BfK+) 






Br{B+- 






Br(B+^B° s K+) ^ 




B s 


Br{B+- 




fpF Bc-*B s 


Br(B+->B° s K*+) 


f K *F B <- 


+B S 


Br(B+- 


-fl°7r+) r 


T7>B C — >B d 
J* 


Br{B+^B° s K + ) 






Br(B+- 


^Bfn+y 


^0 


Br(B+->B?K+) 
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IV. SUMMARY AND CONCLUSION 



In prospects of the huge statistics of the B c mesons at the hadron colliders, accurate and 
thorough studies of the B c physics will be accessible very soon. In this paper, we study the 
two-body nonleptonic c-quark decays in the B c mesons, i.e. B c — > B^P, B q V decays within 
the QCDF approach for the leading approximation, and estimate their branching ratios. We 
find that the contributions of the penguin operators are very small to the decay amplitudes 
due to the serious suppression by the CKM elements. The decay modes determined by a± 
have comparatively large branching ratios. The most promising decay modes are B c — > 
B^n, B s p, which might be easily detected at the hadron colliders. 



APPENDIX A: AMPLITUDES FOR B c -► B { q * } P, B q V DECAYS 



1. c — > d processes 



A(B+^B° d n+) = -^/X^ « - O {K^>i 

1 1 

+ {V,dV c ' d + V U .,V C ;) [«4 - 2«io + " c i(«6 - jag)]} (Al) 



where R cl 



(m d + m u )(m c -m d )' 
A(B^B° d K+) = -i^f K F^ « - m 2 Bd ) V us V^a x (A2) 

A(B+^B° dP +) = V2G F f p F^m p+ (e-p Bc ) {K d K>i 

+ (v ud V; d + V US V; S ) [a 4 - l -a w ] } (A3) 

A(B^B° d K* + ) = V2G F f K *F^ B «m K « + (e-p Bc ) V us V^ ai (A4) 

A(Bt^B* d \ + ) = V2G F f n A B ^ B «°m B * (e-p Bc ) 

+ (v ud V* d + V US V* S ) [04 - ^aio + Q cl (a 6 - ^a 8 )] } (A5) 

Wh6re Qcl = (m d + mu )(m c + md y 

A{Bt^BfK + ) = V2G F f K A B ^ B «°m B * (e-p Bc ) K,K>i (A6) 
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2. c — > s processes 



.G F , pB^B? / 2 



V2" 



m R — m- 



+ + V u S V* s ) [a 4 - ^aio + i^K - ^ 



where i?, 



c2 



2m^+ 



(m s + m u )(m c - m s ) 



A(B+^B° sP + ) = V2G F f p F^ Be m p+ (e-p Bc ) V ud V^ ai 



V2G F f K *F B ^ B "m K , + (e-p Bc ) {v w V* ai 

+(v ud v; d + v us v; s ) [a A - l -a w ]} 

-- V2G F UA^ Bf m Bi (e-p Bc ) V ud V* s a! 



A{Bt^BfK + ) = V2G F f K A B ^ m B% (e-p Bc ) {K,K>i 

+ (V ud V* d + V US V* S ) [a 4 - ia 10 + <2 c2 (a 6 - t^s)] } 



where Q c2 



-2m 2 K+ 



(m s + m u )(m c + m s 



3. c — ► u processes 



(AT) 



} (A8) 



(A9) 



(A10) 
(All) 



(A12) 



A(B^B^) = -i^fUF^ ( 



m R — m' 



? J { - v ud v; d a 2 



where R, 



+ (V ud V* d + V US V*^ [a 4 + aio - -(a 7 - a 9 ) + # c3 (a 6 + a 8 )] } 
2m2 



c3 



(m d + m u )(m c - m u ) 



A(B^B+K°) = -&f K F B ^ (m 2 Bc - m\) V ud V* s a 2 (A14) 



v/2 



(A13) 
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AiB^B^K ) = -i9Zf K F B ^ (m| c - ml u ) V us V c * d a 2 



A(B+^B+ V M) = -i^ (l) F B ^ Bi (m| c - ml u ) {V ud V; d a 2 
+^V us V*a 2 + (V ud V; d + V US V; S ) \^{a 3 -a, + \(a 7 - a 9 )} 

Jr,C) J v (>) A 

+2(a 3 - a 5 ) + a 4 + ai - - (a 7 - a 9 ) + ( 1 - \ i?2(a 6 + a 8 )]} 



where i?c 4 



(m s + m s )(m c - m u ) 

A(B+^B+p°) = G F f p F 1 B ^ Bi m p0 (e-p Bc ){-V vd V^a 2 

+(v ud v^ + v us v c 



3 
2 1 



where Q C 3 = 



A(B+^B+u) = GpUF^m^ (e-p Bc ) {v ud V* d a 2 

+ (V ud V c * d + V US V*) [2(a 3 + a 5 ) + a 4 + a 10 + i(a r + a 9 )] } 

= y/TGFfeFf^mjr* (e-p Bc ) V ud V* s a 2 
A(B+^B+K*°) = V2G F f K ,oF B ^ Bi m K , (e-p Bc ) V us V c * d a 2 

A(B^B: + n°) = G F UA B ^ + m B * (e-p Bc ) { - V ud V c * d a 2 
+ (V ud V* d + V us V*^j [a 4 + aio - -(a 7 - a 9 ) + QdjK + as)] } 

(m d + m u )(m c + m u ) 
^(S+^^+X ) = v / 2G F /^<^ B " + m B:+ ( £ .p flc ) V ud V; s a 2 

A(B^B:+K°) = V2G F f K A B ^ K+ m K+ (e-p B ) V us V c * d a 2 

A{Bt^B* u + ^) = V2G F f^A B ^ + m K+ (e-p B ) {v ud V; d a 2 
+^V us V c * s a 2 + (V ud V? d + V US V; S ) [^{a 3 - a 5 + \ (a 7 - a 9 ) } 



fu 
J v O) 



+2(a 3 - a 5 ) + a 4 + a 10 - - (a 7 - a 9 ) + 1 - Qci (a 6 + a 8 ) } 



where Qc 4 = 



(m s + m s )(m c + m u ) 
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TABLE I: The NLO Wilson coefficient 


3 Cj(/i) in the NDR scheme. The input parameter 


s are 8 ] : 


aJmv) = 1176 


n , ~™(?7?,wl = 1 /1 28 


mw = 80 403 GeV 


A (/=5) 


= 220 9 MeV Al^7^ ) 


= 317.2 


MeV. 














p — '"-6 


,, _ 9 n r; P v 

Lt Zj.\J V_ I c V 


/'' 


— 1 ^ deV 

— J. .tj OC V 






1 D84Q 


1 1AQ7 




1 1 

1 . lOOO 


1.2215 




— D 1 Q09 


— ^077 




— ^71 7 


-0.4241 


(1? 










0.0349 




— n 0369 


—0 0^49 




— 0659 


-0.0747 




n oo&& 






01 07 


0.0102 


c 6 


-0.0422 


-0.0703 




-0.0896 


—0.1078 


Cr/a em 


-0.0007 


-0.0164 




-0.0186 


—0.0181 


Cs/a e m 


0.0565 


0.0964 




0.1235 


0.1493 


Cg/a em 


-1.3039 


-1.3966 




-1.4473 


-1.4901 


Cio/ ct em 


0.2700 


0.4144 




0.4964 


0.5656 



TABLE II: values of the decay constant (in the unit of MeV) 



fn fx 


h 


/. 


fp fui 


Ik* 


131 [8] 160 [8] 


(1.07±0.02)/ 7r [44] 


(1.34±0.06)/ 7r [44] 


205±9 [53] 195±3 [53] 


217±5 [53] 



a The form factors increase with the increasing parameter u — 0.4 ~ 1.0 GeV that determines the aver- 
age transverse quark momentum. The authors of [37| prefer F Bc ^ Bu (0) = 0.831, F B ^ B > (0) = 0.859, 
A Bc ^ B: (0) = 0.869 and A Bc ^ B ' s (0) = 0.842 with the corresponding parameter uj = 0.8 GeV. 

b The definitions of the transition form factors in [46j are different from ours in Eq. (]2"4")) and Eg . (f2"S")) . The 
relationship is 

Ff^.F + . (39) 

with the values of F + = 0.3±0.1 (0.30±0.05), = 4.0±1.0 (4.5±0.5) GeV^and = -0.02±0.01 

(-0.03±0.02) GeV" 1 for B c -> B ( * ] d (B { s * ] ) transition [|6|. 
c Using the relationship of Eq.® and Eq.flHO} with the input A x = 0.52, A 2 = -2.79 [33j. 
d For parameter ui = 0.4, 0.5 GeV. 

c Using the relationship of Eq.® with the input F + = 0.4504 (0.5917), F A = 3.383 (5.506) GeV" 1 and 

F A = -0.0463 (-0.0673) GeV" 1 for B c B^ d (B^) transition Q- 
f Using the relationship of Eq.® and Eq.® with the input A x = 0.27 (0.33) and A 2 = -0.60 (-0.40) 

for B c -> B { *\ (B^) transition Q- 
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TABLE III: Values of transition form factors 



Ref. 


F Bc " B,1 ' d (0) 




B r — *B i , 

A, - d (o) 


A*°-* B *(0) 


[37] a 


0.320~0.910 


0.340~0.925 


0.349~0.916 


0.432^0.931 


[46] h 


0.3±0.1 


0.30±0.05 


0.35±0.09 


0.39±0.05 


[33J c 




0.61 




0.79 


[47] d 




0.403~0.617 




0.433~0.641 


[48] e 


0.4504 


0.5917 


0.2691 


0.4451 


[49] / 


-0.58 


-0.61 


0.35 


0.39 


[50] a 




0.297 




0.263 


[21] fc 


1.27 


1.3 


1.29 


0.94 


[21] i 


1.38 


1.1 


1.26 


1.04 


[17J 


0.39 


0.50 


0.20 


0.35 


[5JJ 






0.23±0.03 




[52] ; 


0.90 


1.02 


0.27 


0.36 



s Using the relationship of Eq.{29]) and Eq.(30]) with the input A x = 0.28 and A 2 = 0.49 0j- 

h Using the relationship of Eq.® with the input F+ = 1.27 (1.3), F A = 9.8 (8.1) GeV -1 and = 0.35 

(0.2) GeV -1 for B c -> B { *\ (B^) transition in the framework of QCD sum rules [2l| . 
1 Using the relationship of Eq.® with the input F + = 1.38 (1.1), F A = 9.4 (8.2) GeV^ 1 and Fjj 4 = 0.36 

(0.3) GeV -1 for B c — ► B^*\ (B^) transition in the framework of potential model [21] . 
j Using the relationship of EgM)§ and Eq.j30|) with the input A 1 = 0.90 (1.01) and A 2 = 7.9 (9.04) for B c 

-> B { * ] d {Bi* ] ) transition Q. 
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TABLE IV: The branching ratios for B c — > P, B q V . Br T corresponds to the contributions 
of the operators Q\ and Q2- Br T+Ps corresponds to the contributions of operators Q\ ~ Q§. 
q t t+p 3 +p c correS ponds to the contributions of Q\ ~ Qio- 



modes case Br T Br T + p ° Br T + p *+ p * BrT+ ?,f BrT Br^+^-Br? 





~* 




case la 


5.3089xl(T 


-2 


















B+ 


~* 


B° sP + 


case la 


6.2652x10" 


-2 


















Bt 


-* 


Bfir+ 


case la 


4.5916x10- 


-2 


















Bt 


-* 


B° S K+ 


case lb 


3.6746x10- 


-3 


3.6759x10" 


-3 


3.6759x10" 


-3 


3.4x10" 


-4 


3.4x10" 


-4 


Bt 


-* 


B° S K*+ 


case lb 


1.6450x10- 


-3 


1.6451x10" 


-3 


1.6451x10" 


-3 


5.0x10" 


-5 


5.0x10" 


-5 


Bt 


-* 


BfK+ 


case lb 


2.9772x10- 


-3 


2.9766x10" 


-3 


2.9766x10" 


-3 


-1.9x10" 


-4 


-1.9x10" 


-4 


Bt 


~* 


B> + 


case lb 


3.7283x10- 


-3 


3.7272x10" 


-3 


3.7272x10" 


-3 


-3.0x10" 


-4 


-3.0x10" 


-4 


Bt 


-* 


B° d P + 


case lb 


5.2745x10- 


-3 


5.2742x10" 


-3 


5.2742x10" 


-3 


-5.0x10" 


-5 


-5.0x10" 


-5 


Bt 


-* 


Bfn+ 


case lb 


3.2682x10- 


-3 


3.2688x10" 


-3 


3.2688x10" 


-3 


1.9x10" 


-4 


1.9x10" 


-4 


Bt 


-* 


B° d K+ 


case lc 


2.6616x10- 


-4 


















Bt 


-* 


B° d K* + 


case lc 


2.2583x10- 


-4 


















Bt 


-* 


B* d °K+ 


case lc 


2.2075x10- 


-4 


















Bt 


~* 


B+K° 


case 2a 


2.2067x10- 


-5 


















Bt 


-* 


B+K*° 


case 2a 


1.8434x10' 


-5 


















Bt 


~* 


B*+IC 


case 2a 


1.8261x10" 


-5 


















Bt 


-* 


B+V 


case 2b 


1.5991x10- 


-6 


1.6122x10" 


-6 


1.6125x10" 


-6 


8.2x10" 


-3 


8.4x10" 


-3 


Bt 




B* + V 


case 2b 


1.3042x10- 


-6 


1.2960x10" 


-6 


1.2964x10" 


-6 


-6.3x10" 


-3 


-6.0x10" 


-3 


Bt 




BU° 


case 2b 


4.5968x10- 


-7 


4.5161x10" 


-7 


4.5134x10" 


-7 


-1.8x10" 


-2 


-1.8x10" 


-2 


Bt 




B+P° 


case 2b 


6.5030x10' 


-7 


6.4823x10" 


-7 


6.4776x10" 


-7 


-3.2x10" 


-3 


-3.9x10" 


-3 


Bt 




B+u 


case 2b 


5.7921x10- 


-7 


5.8199x10" 


-7 


5.8212x10" 


-7 


4.8x10" 


-3 


5.0x10" 


-3 


Bt 






case 2b 


4.0262x10- 


-7 


4.0722x10" 


-7 


4.0685x10" 


-7 


1.1x10" 


-2 


1.0x10" 


-2 


Bt 




sir/ 


case 2d 


8.8676x10- 


-8 


8.7700x10" 


-8 


8.7738x10" 


-8 


-1.1x10" 


-2 


-1.1x10" 


-2 


Bt 




B* + V' 


case 2d 


1.7401x10- 


-8 


1.7728x10" 


-8 


1.7731x10" 


-8 


1.9x10" 


-2 


1.9x10" 


-2 


Bt 




B+K° 


case 2c 


6.5428x10- 


-8 


















Bt 




B+K*° 


case 2c 


5.4658x10- 


-8 


















Bt 




B*+K° 


case 2c 


5.4143x10- 


-8 
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